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ABSTRACT 

The young composite supernova remnant (SNR) Kesteven 75, with a pulsar wind nebula at its 
center, has an unusual morphology with a bright southern half-shell structure in multiwavelengths. 
The distance to Kes 75 has long been uncertain. Aiming to address these issues, we have made 
millimeter spectroscopic observations of the molecular gas toward the remnant. The Vlsr ~ 83- 
96 km s _1 molecular clouds (MCs) are found to overlap a large north-western region of the remnant 
and are suggested to be located in front of the SNR along the line of sight. Also in the remnant 
area, the Vlsr= 45-58 km s _1 MC shows a blue-shifted broadening in the 12 CO (J=l-0) line profile 
and a perturbed position- velocity structure near the edge of the remnant, with the intensity centroid 
sitting in the northern area of the remnant. In particular, a cavity surrounded by a molecular shell is 
unveiled in the intensity map in the broadened blue wing (45-51 km s _1 ), and the southern molecular 
shell follows the bright partial SNR shell seen in X-rays, mid- infrared, and radio continuum. These 
observational features provide effective evidences for the association of Kes 75 with the adjacent 
54 km s" 1 MC. This association leads to a determination of the kinematic distance at ~ 10.6 kpc to 
the remnant, which agrees with a location at the far side of the Sagittarius arm. The morphological 
coincidence of the shell seen in multiwavelengths is consistent with a scenario in which the SNR shock 
hits a pre-existing dense shell. This dense molecular shell is suggested to likely represent the debris 
of the cooled, clumpy shell of the progenitor's wind bubble proximately behind the 54 km s _1 cloud. 
The discovery of the association with MC provides a possible explanation for the 7-ray excess of the 
remnant. 

Subject headings: ISM: individual (Kes 75, G29.7— 0.3) - ISM: molecules - supernova remnants 



1. INTRODUCTION 

Massive stars have often not moved far from their ma- 
trices in molecular clouds (MCs) by the time they ex- 
plode as core-collapse supernovae because of their short 
lifetime, and therefore their remnants are often observed 
to be associated with MCs. The CO observations have 
demonstrated a good positional correlation between MC 
complexes and supernova remnants (SNRs) (Huang & 
Thaddeus 1986). About 20 SNRs are known to phys- 
ically interact with ambient MCs with convincing evi- 
dences (e.g. Frail et al. 1996). The interaction of either 
the progenitors' stellar winds or the SNR shocks with 
the MCs can play an important role in the SNRs' evo- 
lution and morphologies. The obscuration by the MCs 
can also heavily affect the X-ray morphologies. Some of 
SNRs have been known to have a peculiarly asymmetric 
morphology in multiwavelengths because of the associ- 
ation with MCs. For example, the western blast wave 
of SNR CTB 109 is thought to have been significantly 
decelerated by a dense MC, resulting in an eastern semi- 
circular shape of the remnant (Tatematsu et al. 1987; 
Wang et al. 1992; Rho & Petre 1997; Sasaki et al. 2004). 
The soft X-rays of the northwestern part of SNR 3C 391 
are substantially extincted by the MC in which the SNR 
is embedded (Rho & Peter 1996; Wilner et al. 1998; 
Chen & Slane 2001; Chen et al. 2004). Following this 
notion, other SNRs with highly asymmetric (especially, 
half-brightened) morphologies, such as Kesteven 75 and 
Kesteven 69 (Zhou et al. 2008; hereafter Paper I), are 



thus of particular interests. 

Kesteven 75 (G29.7— 0.3) is one of the composite SNRs, 
exhibiting both the Crab-like and shell-type properties 
(Becker & Helfand 1984; Blanton & Helfand 1996). A 
young energetic X-ray pulsar, PSR J1846-0258, discov- 
ered using the Rossi X-ray Timing Explorer (RXTE), 
was localized to within an arcminute of the remnant cen- 
ter using an ASCA observation (Gotthelf et al. 2000). 
The bright point-like source at the core of the remnant 
was later resolved with the high resolution Chandra ob- 
servation (Helfand et al. 2003). Notably, both the Very 
Large Array (VLA) 1.4 GHz radio and Spitzer 24/im mid- 
infrared (IR) observations show a semicircular partial 
shell (3.5' in diameter) in the south (Becker & Helfand 
1984; Morton et al. 2007). The Chandra X-ray image 
of Kes 75 clearly exhibits a synchrotron-emitting pulsar 
wind nebula (PWN) in the center, two bright elongated 
patches (in the southeast and southwest, respectively) 
along the radio and IR partial shell, and very faint dif- 
fuse emission in the northern part (Helfand et al. 2003). 
Recently, by deep Chandra X-ray observations, the X-ray 
softening and temporal flux enhancement of the central 
pulsar were revealed (Kumar & Safi-Harb 2008), and a 
detailed spatial and spectral analysis of the PWN was 
made (Ng et al. 2008); in both works, a distance 6 kpc 
was used. 

On the other hand, the distance to Kes 75 is under 
controversy. By HI absorption observations, Caswell 
et al. (1975) gave a range ~ 6.6-19 kpc for the kine- 
matic distance, covering the value ~ 19 kpc derived by 
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Milne (1979) from the £-D relationship, while Becker 
& Helfand (1984) suggested - 21 kpc. With a higher- 
sensitivity HI observation, Leahy & Tian (2008, here- 
after LT08) provided a small estimate of ~ 5.1-7.5 kpc. 
McBridc ct al. (2008) favor a distance as small as ~ 3 kpc, 
so as to reconcile the problem of the otherwise very large 
efficiency in converting the central pulsar's spin-down 
power into X- and 7-ray luminosity. 

Investigating the physical relation between Kes 75 and 
its ambient MCs may be very helpful to understanding 
its evolution and irregular morphology as well as clarify- 
ing its distance. In this paper, we present our millimeter 
CO observations toward this remnant. In § 2, we describe 
the observations and the data reduction. The main ob- 
servational results are presented in § 3, and discussion 
and conclusions are in § 4 and § 5, respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

The observation was first made in the 12 CO (J=l-0) 
line (at 115.271 GHz) in 2006 April using the 6-meter 
millimeter-wavelength telescope of the Seoul Radio As- 
tronomy Observatory (SRAO) with single-side band fil- 
ter and the 1024-channel autocorrelator with 50-MHz 
bandwidth. The half-power beam size at 115 GHz is 
2' and the typical rms noise level was about 0.2 K at 
the 0.5 km s _1 velocity resolution. The radiation tem- 
perature is determined by Tr = Tx/(fb x Vmb), where 
Ta is the antenna temperature, fb the beam filling fac- 
tor (assuming ft ~ 1), and rj^ the main beam efficiency 
(~ 75%). We mapped the 10' x 10' area covering Kes 75 
centered at (18 h 46 m 25 s .0, -02°59'00") with grid spacing 
~ T. 

The follow-up observations were made in the 
12 CO (J=l-0) line (at 115.271 GHz), 13 CO (J=l-0) line 
(110.201 GHz), and C 18 (J =1-0) line (109.782 GHz) 
during 2006 November-December and 2007 October- 
November using the 13.7-meter millimeter-wavelength 
telescope of the Purple Mountain Observatory at Del- 
ingha (hereafter PMOD). The spectrometer has 1024 
channels, with a total bandwidth of 145 MHz (43 MHz) 
and the velocity resolution of 0.5 km s _1 (0.2 km s _1 ) 
for 12 CO ( 13 CO and C 18 0). The half-power beamwidth 
of the telescope is 54" and the main beam efficiency 
Vmb is 67%. We mapped the source centered at 
(18 h 46 m 25 s .0,-02°59'18") with grid spacing 0.5' in the 
inner (4.5' x 5') region (rms noise < 0.1 K) and 1' in 
the outer region (out to 10' x 10', rms noise ~ 0.3 K). 
All the CO observation data were reduced using the 
GILDAS/CLASS package developed by IRAM 1 . 

The Chandra X-ray and Spitzer 24 ^m mid-IR obser- 
vations are also used. We revisited the Chandra ACIS 
observational data of SNR Kes 75 [ObsIDs 748 (PI: D.J. 
Helfand), 6686, 7337, 7338, & 7339 (PI: P. Slanc) with 
total exposure time of 196 ks]. We reprocessed the event 
files (from Level 1 to Level 2) using the CIA03.4 data 
processing software to remove pixel randomization and 
to correct for CCD charge-transfer inefficiencies. The 
overall light curve was examined for possible contamina- 
tion from a time- variable background. The reduced data, 
with a net exposure of 183 ks, were used for our analysis. 
The mid-IR 24 fim observation used here was carried 
out as the program of the Breaking Down the Spectra 

1 http: / /www.iram.fr /IRAMFR / GILDAS 



of Pulsar Wind Nebulae (PID: 3647, PI: P. Slane) with 
the Multiband Imaging Photometer (MIPS) (Rieke et al. 
2004). The Post Basic Calibrated Data (PBCD) of the 
24 /zm mid-IR are obtained from Spitzer archive. The 
1420 MHz continuum and Hl-line data come from the 
archival VLA Galactic Plane Survey (VGPS) (Stil et al. 
2006). 

3. RESULTS 

3.1. Velocity Structure of Molecular Gas Components 
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Fig. 1. — 12 CO (J= 1-0) spectrum from the SRAO observation 
(0 < V < 120 km s" 1 ) toward SNR Kes 75 in a 10' X 10' area 
centered at (18 h 46 m 25 s .O, -02°59'00"). In the little box is the 
13 CO(J=1-0) spectrum from the PMOD observation (20 < V < 
120 km s- 1 ) toward SNR Kes 75 in a 4.5' X 5' area centered at 
(18 h 46 m 25 s .0, -02°59'18"). 

We show in FigureQ]the 12 CO (J=l-0) spectrum from 
the SRAO observation in the velocity range 0-120 km s _1 
toward Kes 75 in a 10' x 10' area. Beyond this range, 
virtually no 12 CO emission is detected at —110 km s _1 < 
Vlsr < km s _1 and Vlsr > 115 km s _1 in the 
PMOD observation. The CO emission in the direction of 
SNR Kes 75 is present in a broad velocity range between 
3 and 110 km s _1 , characterized by several prominent 
peaks. The first 12 CO component is represented by the 
primary emission peak in the interval Vlsr=3-17 km s _1 . 
We call the molecular gas at 40-60 km s _1 as the sec- 
ond component and will show its association with the 
SNR. The other 12 CO emission peaks in the interval 60- 
110 km s _1 are referred to as the third component of 
MC complex, although the peaks in this interval are not 
necessarily related to each other. The first and the third 
velocity components are probably corresponding to the 
Aquila Rift (Dame et al. 2001) and the MCs in the Scu- 
tum and 4-kpc arms (Dame et al. 1986), respectively. 
The second velocity component displays a narrow Gaus- 
sian profile interval of 49-58 km s _1 around ^54 km s _1 
and an additional broad line emission structure interval 
of 45-51 km s -1 (with an average 6 km s -1 blueshift). In 
the 13 CO (J=l-0) observation (20 < V < 120 km s" 1 ), 
there are the counterparts of the second and third molec- 
ular components. No signal is detected in the range 20- 
120 km s" 1 in the PMOD C 18 ( J=l-0) observation. 

We made the 12 CO-emission intensity channel maps 
with ~ 1 km s _1 velocity intervals. No evidence is found 
for the positional correlation between the first CO com- 
ponent and SNR Kes 75. However, the second and third 
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TABLE 1 

Parameters of the Gaussian profile at 
Vlsr=49-58 km s _1 derived from the PMOD 
observation. 



18 h 46 m 42 s 36 s 30 s 24 s 18 s 12 s 18 h 46 m 42 s 36 s 30 s 24 s 18 s 12 s 18 h 46 m 42 s 36 s 30 s 24 s 18 s 12 s 

RA (J2000) RA (J2000) RA (J2000) 

Fig. 2.— The SRAO 12 CO (J=l-0) intensity channel maps between 49 km s 1 and 58 km s 1 (smoothed to a resolution of 0.2' by 
interpolation), overlaid with the 1.0-7.0 keV X-ray (cyan, thin) and 1.4 GHz radio (red, thick) continuum emission contours (from the 
NRAO VLA Sky Survey). 

components seem to have positional correlations with 
the remnant. Figures [2] and [3] show the 12 CO channel 
maps in the intervals 49-58 km s _1 and 76-101 km s _1 
(the main bodies of the second and third velocity compo- 
nents), respectively, with 1' grid pointing. The channel 
maps are overlaid with the contours of the Chandra 1.0- 
7.0 keV X-ray and VLA 1.4 GHz radio continuum emis- 
sion, respectively. In morphology, both the X-ray and 
radio emissions show bright half in the south and faint 
half in the north. It is interesting that the 49-58 km s _1 
cloud (as a part of the second component) is projection- 
ally situated north to the radio- and X-ray-bright south- 
ern half of the SNR, seeming to cover the faint northern 
half. The CO emission around 51 km s _1 also appears to 
surround the SNR from the north to the west. As parts 
of the third component of MC complex, the Vlsr ~ 83- 
96 km s" 1 MCs overlap a large north-western region of 
Kes 75, and the ~ 97-101 km s _1 MCs completely over- 
lap the remnant region (Fig. [3]). The brightness anti- 
correlation between the two 12 CO line components and 
the X-ray and radio (as well as mid-IR 24^m, §3.2|) emis- 
sions is clear in a large scale. 
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a The half velocity width of the line. b The region of the MC 
in the interval 49-58 km s _1 is about 27 arcmin 2 . 

3.2. The Vlsr ~ 54 km s _1 Molecular Cloud 
Component 

In order to obtain detailed clues to the possible associa- 
tion of the cloud components with the SNR, we have ana- 
lyzed the second and third velocity components in the in- 
ner 4.5' x 5' area [centered at (18 h 46 m 25 s .O, -02°59'18")] 
based on the PMOD observation. 

The observed parameters of the CO lines and the de- 
rived column density and cloud mass for the second ve- 
locity component (in interval 49-58 km s -1 ) are summa- 
rized in Table 1 and Table 2, respectively. Two meth- 
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Fig. 3.— The SRAO 12 CO(J=1-0) intensity channel maps between 76 km s 1 and 101 km s 1 (smoothed to a resolution of 0.2' by 
interpolation), overlaid with the 1.0—7.0 keV X-ray (cyan, thin) and 1.4 GHz radio (red, thick) continuum emission contours (from the 
NRAO VLA Sky Survey). 



ods have been used in the derivation of the molecu- 
lar column density and similar results are yielded. In 
the first method, the H2 column density is calculated 
by adopting the mean CO-to-H2 mass conversion fac- 
tor 1.8 x 10 20 cm _2 K -1 km _1 s (Dame et al. 2001; sim- 
ilar values also derived by Strong & Mattox 1996 and 
Hunter et al. 1997). In the second method, on the as- 
sumption of local thermodynamic equilibrium (LTE) and 
12 CO (J=l-0) line being optically thick, the 13 CO col- 
umn density is converted to the H2 column density us- 
ing 7V(H 2 )/iV( 13 CO) w 7 x 10 5 (Frerking et al. 1982). 
In the estimate of the LTE mass of the molecular gas, 
a mean molecular weight per H2 molecule 2.76toh has 
been adopted. Assuming the line of sight (LOS) size of 
the cloud similar to its apparent size ~ 4', we estimate 



the molecular density as n(H2) 



60d 



1 — s 
cm 3 



10.6 



where 



d W 6 = d/(10.6 kpc) is the distance to the MC/SNR in 
units of the reference value (see H4.ll for the distance de- 
termination) . 

We extracted the 12 CO and 13 CO spectra in the ve- 



locity range 40-60 km s _1 from two selected regions 
(as defined in Fig. [4]) , which basically correspond to the 
X-ray/IR/radio-bright southern half and faint northern 
half, respectively. The 12 CO and 13 CO spectra (shown 
in Fig. Ej) in the both regions display a Gaussian profile 
around Vlsr ^54 km s _1 . There is a broad left wing (45- 
51 km s _1 ) in the 12 CO line profiles for both the northern 
and southern regions (as already seen in the spectrum for 
the entire remnant field, Fig. [1]). Similar wing broaden- 
ing is not seen in the 13 CO lines. The ~ 49-58 km s _1 
12 CO and 13 CO Gaussian lines of the northern region 
are stronger than those of the southern region, indicating 
that the CO emission come majorly from the northern 
region, consistent with the 12 CO channel maps (Fig. [2]). 
However, the broadened left wing is not as strong in the 
northern region as in the southern one except in the in- 
terval 49-51 km s _1 . For the defined southern region, 
we obtain the mean integrated intensity spanning the 
line profile (45-58 km s _1 ) ~8.2 K km s _1 and the cor- 
responding molecular column density and gas mass as 
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TABLE 2 

Parameters for. the MC derived from the PMOD observation. 



Region Vlsr interval 
(km s _1 ) 


H2 column density a 
N(H 2 ) (10 2Q cm- 2 ) 


Virial mass b 
M vir (1O 3 M ) 


Molecular mass 
M (1O 3 M ) 


Main body 49-58 


22.4 / 21.5 


26.6dio. 6 c 


12.6d? n „ / 12.1d? n fi ° 


Southern region d 45—58 
Gaussian around 54 km s — 1 
Gaussian subtracted 


14.8 

9.6 

5.2 




2-7^ .6 C c 

1 - 8d 10.6 C 

0.9d? n K c 



a See text for the two methods used for the column density derivation. 11 The virial mass is obtained 
using the formula M vir = 210(r/lpc)(Ai;/lkm s -1 ) 2 Mq, where Av is the FWHM of the line 
(Caselli et al. 2002). c Parameter dio.a is the distance to the cloud in units of 10.6 kpc ( i|4.1> . d A 
sky region about 9 arcmin 2 , see Fig|4] 
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Fig. 4. — The PMOD 12 CO (J=l-0) intensity maps in the intervals 45-47, 47-49, 49-51, & 45-51 km s -1 . The two regions (6.5 arcmin 2 
upper and 9 arcmin 2 lower) in the upper-right panel are defined for the CO-spectrum extraction. In the lower-right panel, the 12 CO 
intensity map is overlaid with the 6 cm radio continuum contours (Becker & Helfand 1984), the star stands for the position of the pulsar, 
and the dotted lines indicate the E-W and N-S cuts for the PV diagrams (Fig. [6jl . 



listed in Table 2. The contribution of the Gaussian com- 
ponent centered at 54 km s _1 and the broadened part 
(with the Gaussian component subtracted) are also esti- 
mated, respectively, as given in Table 2. 

To better understand the velocity structure of the 
second 12 CO component, we made the position- 
velocity (PV) diagrams (Fig. [6]) along the east- 
west (E-W) and north-south (N-S) cuts centered at 



(18 h 46 m 25 s .0,-02°58'48") (as labeled in Fig. g]). In Fig- 
ure [6j it is clearly shown that the main body of the MC 
is located in the north of the Kes 75 area and the maxi- 
mums of the 12 CO line broadening appear near the edge 
(at radius about 1.7-2' from the pulsar) of the remnant. 

The close-up 12 CO channel maps and the integrated 
map in the blue- wing interval 45-51 km s _1 are presented 
in Figure 21 in which the 6-cm radio continuum contours 
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Fig. 



5. — The 40—60 km s CO spectra from northern 
(6.5 arcmin 2 , left) and southern (9 arcmin 2 , right ) regions of SNR 
Kes 75 based on the PMOD observation. The thick and thin lines 
are for the 12 CO(J=1-0) and 13 CO(J=1-0) emission, respectively. 

are superposed on the integrated map. In Figure [7] is 
shown a tri-color image of the Kes 75 area, with the 
12 CO intensity in the interval 45-51 km s _1 in red, the 
Spitzer 24{im mid-IR intensity (together with the con- 
tours) in green, and the Chandra X-ray emission in blue. 
A cavity is revealed to be surrounded by a molecular 
shell in the line wing (45-51 km s _1 ) intensity map, and, 
surprisingly, the shell in the south is coincident with the 
bright rim seen in the X-ray, mid-IR, and radio emis- 
sions. This shell is responsible for the maximums of the 
line broadening in the PV diagrams. The bright 12 CO 
patch on the northern edge of the SNR may contain the 
contamination from the contribution of the main body of 
the ~ 54 km s _1 MC at its Gaussian wing 49-51 km s _1 
(left panel of FigEJ) , which has been seen to be located 
in the north (Fig[2]and FigJS]). 

We did not find any CO feature between 60 km s _1 
and 110 km s _1 (i.e. for the third molecular component) 
positionally coincident with the southern X-ray/mid- 
IR/radio-bright rim. The clouds in this velocity range 
are estimated to have a column density iV(H2) ~ 2.4 x 
10 22 cm~ 2 in the northern region (as defined in Fig[4| 
and ~ 1.9 x 10 22 cm~ 2 in the southern region (Figf5|). 



3.3. X-ray Spectral Properties of Kes 75 

To investigate the positional relation of the X-ray prop- 
erties of SNR Kes 75 with the CO intensity distribution, 
we extracted three X-ray spectra of the southern clumpy 
shell and the faint northern diffuse emission (Fig. [8]) 
from the reduced Chandra data using the CIA03.4 script 
specextract. The background spectrum was extracted 
from a southeastern region outside the shell. The three 
on-source spectra are adaptively regrouped to achieve a 
background-subtracted signal-to-noise ratio (S/N) of 5 
per bin. The previous Chandra X-ray studies by Helfand 
et al. (2003) and Morton ct al. (2007) both used a uni- 
form absorbing hydrogen column density Nh, fixed to 
that for the central PWN, for the entire remnant. Here 
we will allow Nh to vary for each region in our spectral 
fit, in which the XSPEC11.3.2 software package with the 
Morrison & McCammon (1983) interstellar absorption is 
used. 

The spectra of the two clumps along the southern shell 
contain distinct Hea lines of Mg, Si, S, and Ar (see FigJ5] 
for the southeastern clump), as have been seen in the 
previous studies. However, the spectrum of the northern 
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Fig. 6. — Position-velocity diagrams along E-W cut (upper panel) 
and N-S cut (lower panel) [centered at (18 h 46 m 25 s .0, -02°58'48")] 
based on the PMOD observation. Contour levels start from 1.0 K 
with in steps of 0.4 K. The dashed lines mark spatial the regi on o f 
radius 1.7' from the pulsar and the thick line marks 2.0' (see i|4,2l l. 
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Fig. 7.— Tri-color image of SNR Kes 75. The PMOD 12 CO 
intensity map (> 2.2 K km s — 1 or llcr) in the interval 45—51 km s — 1 
is in red, the Spitzer 24/xm mid-IR emission in green (including 
green contours), and the Chandra 1.0—7.0 keV X-rays (adaptively 
smoothed with the CIAO program csmooth to achieve a S /N ratio 
of 3) in blue. Each of the intensity maps is in logarithmic scales. 



region appears to be featureless (Fig [TU|) . The northern 
diffuse emission may be the dust-scattered light from the 
PWN; actually, Helfand et al. (2003) have pointed out 
that the effect of the dust-scattering halo of the Kes 75 
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TABLE 3 

Chandra X-ray Spectral Fits with 90% Confidence 
Ranges 



Parameter 


SE Clump 


SW Clump 


Northern Diffuse 


A ff (10 22 cm" 2 ) 


q p.Q + 0.05 
°- u,:> -0.13 


Q '7Q _ j _ 0.20 

z * ' z -0.17 


Q Q1 +0.47 

-0.32 


kT x (keV) 


n ^r + 03 

- 65 Iao2 


°- 88 -0 21 




Mg 


1 „„+0 18 






Si 


i 7 q+0.28 
1 -' a -0.26 


ASS -0.83 




S 


9 7C+0.37 
z - '°-0.40 


7 lfi+ 4 ' 57 




r(10 11 s cm" 3 ) 


4 12+ 1 ' 23 


f) QC+0.10 

U - JO -0.07 




norm(l0~ 2 ) 


i oc+0.35 
1,JU -0.44 






Pholndex 


9 7R+0-29 


2.801^7 


9 7^+0.29 
z - ,o -0.20 


norm (10 ) 


l.92t°;^ 


1 54+O.I8 
1-0 -0.14 


n c fi +0.28 
U ' OD -0.14 


X 2 /dof 


1.19/277 


1.14/244 


0.86/81 


A ff (10 22 cm" 2 ) 


n r,n + U.U2 

°'^°-0.05 


3.55_ n 




kT x (keV) 


n en + 0.04 
°- 60 -0.03 


S2+ ' 21 
u - oz -0.06 




r(10 12 s cm" 3 ) 


> 2.8 


> 1.0 




norm(10" 2 ) 


1 84+ ' 28 


3.97± 2 '^ 




kT x (keV) 


9 09+O.33 


9 91+0:26 

^• zl -0.20 




Mg 


1 Sfi+ ' 27 
i.oo_ _ 26 


< 0.42 




Si 


, fi7 +0.128 
1 - D '-0.16 


1 11+O.I8 
1 - 11 -0.16 




S 


9 oq+0.23 
Z - Os, -0.34 


1 7 9+0.23 
1 -' z -0.25 




r(10 1(, s cm" 3 ) 


t-o+0.36 
°- JO -0 .44 


3 41+ - 56 

-0.37 




norm(10" 3 ) 


9 70+0.47 
2-78_ 046 


oc-+0.57 
■3-JO_o.58 




X 2 /dof 


1.15/277 


1.26/244 





-2:57:00.0 
30.0 
c 58:00.0 
■| 30.0 
£ 59:00.0 
$ 30.0 
°-3:00:00.0 
30.0 
01:00.0 
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Fig. 8. — The Chandra X-ray raw image of SNR Kes 75, labeled 
with the defined regions for the X-ray spectrum extraction. The 
point sources (except the central pulsar) detected with a wavelet 
source-detection algorithm have been removed. 

PWN should be taken into account. Recently, Seward 
et al. (2006) clearly detected a dust-scattered X-ray halo 
around the Crab comprising 5% of the total strength, 
with surface brightness measured out to a radial dis- 
tance of 18'. Therefore we fit the featureless spectrum of 
the northern region with a power-law model and, follow- 
ing Helfand et al. (2003), fit the spectra of the southern 
clumpy shell with a vnei + power-law model. As an al- 
ternative, we also fit the two spectra of the southern shell 
with a double-thermal-componcnt model (accounting for 
the forward- and reverse-shocked gas), as Morton et al. 
(2007) did. The parameters of the fitting (with the 90% 
confidence ranges) are listed in Tabic 3. 

The both models of our spectral fit generate a lower 
intervening hydrogen column density for the southern 
regions (< 3.6 x 10 22 cm" 2 ) than that used in the 




channel energy (keV) 

Fig. 9. — The Chandra-ACIS X-ray spectrum of the southeast- 
ern clump of Kes 75 (taken from a 183 ks data), fitted with the 
vnei (dashed line) plus power-law (dashed-dotted line) model. The 
dotted line around 3.12 keV is a Gaussian accounting for the Ar 
XVII line. 
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Fig. 10. — The Chandra- ACIS X-ray spectrum of the northern 
region of Kes 75 (taken from a 183 ks data), fitted with a non- 
thermal power-law model. 

previous studies (fixed to the value for the PWN ~ 
4.0 x 10 22 cm" 2 ). The northern diffuse emission has an 
intervening hydrogen column Nn ~ 3.5-4.3 x 10 22 cm" 2 , 
which seems higher than that for the southern re- 
gions. The column estimate for the southern regions 
from the vnei+powerlaw model, 2.6-3.1 x 10 22 cm" 2 , 
is clearly lower than the column for the northern re- 
gion. It is notable that the difference of Nh between 
the northern region of the remnant (including the cen- 
tral PWN) and the southern regions in the case of us- 
ing the vnei+powerlaw model is around ~ 1 x 10 22 cm" 2 
(the differnce between the northern and southwestern re- 



gions is l.ljlg 4 x 10 22 cm 2 ancl that between the north- 
ern and southeastern regions is 0.8^04 x 10 22 cm -2 ). 
The both models produce elevated Si and S abundances. 
Using the volume emission measure, fn e riHV, derived 
from the normalization of the XSPEC model, we ob- 
tain an estimate for the hot-gas hydrogen-atom den- 
sity n H ~ 5.4/1.4/-V 2 d 

10.6 cm f° r the southeast- 
ern/southwestern regions, where n e is the electron den- 
sity, V is the volume of the emitting region, and / is 
the volume filling factor of the hot gas (0 < / < 1). 
Here we have assumed n e ss 1.2n# and oblate spheroids 
(24" x 42" x 42" and 19" x 57" x 57", respectively) for the 
defined southeastern and southwestern elliptical regions 
(Fig.[8|). The mass of the southeastern/southwestern X- 

ray bright region is ~ 3.4/1. 7 f x / 2 d^ 6 MQ. The ioniza- 
tion timescales of the southeastern/southwestern X-ray- 
brighten regions in the vnei+powerlaw model thus im- 
ply ages of ~ 2.0/0.6/ 1 / 2 <i 1 Q 2 6 kyr, which are compara- 
ble with the spin-down timescale 723 yr (Gotthelf et al. 
2000), considering factor / is uncertain. 

4. DISCUSSION 

4.1. Association of SNR Kes 75 with the 
VhSR km s" 1 MC and the Distance 

From the above results of CO observations, we have ob- 
tained evidences for the association between SNR Kes 75 
and the Vlsr ^54 km s _1 MC. First, as seen in the 
12 CO spectra (Figs. [T]&[5]), the broadened blue wing of 
the second (54 km s _1 ) MC component indicates that 
the MC contains a part of gas that suffers a pertur- 
bation toward the observer. Similar one-sided broad- 
ened velocity-profiles are seen from a piece of MC in 
SNR CTB 109 (Fig. 2 in Sasaki et al. 2006) and the 
molecular arcs in SNR Kes 69 (Paper I). Secondly, the 
PV diagrams (Fig. [6]) suggest that this blue- wing broad- 
ening takes place chiefly near the edge of the Kes 75 
remnant. These characteristics hint that the 54 km s _1 
MC component may be perturbed by some interaction 
related with SNR Kes 75. Thirdly, the image of this 
broadening part is characterized by a shell-like structure 
surrounding a cavity, the southern part of which is coin- 
cident with the partial SNR shell seen in X-ray, mid-IR, 
and radio continuum (Figs. [4] and [7]). This coincidence 
not only, again, demonstrates a consistent spatial po- 
sition with that of the SNR, but also implies a shock 
interaction between the SNR and the molecular material 
(detailed in p~2]) . 

The establishment of the association of the 
Vlsr ~54 km s^ 1 MC with SNR Kes 75 can provide an 
independent estimate for the kinematic distance to the 
SNR. Corresponding to the line center Vlsr ~54 km s _1 
of the unperturbed Gaussian profile f iJ3.2p . there are 
two candidate kinematic distances to the MC along the 
LOS in the direction of Galactic longitude I = 29.7°. 
Using the rotation curve in Clements (1985), the two 
candidate distances are 3.3 and 10.6 kpc, which fall in 
the near and far side of the Sagittarius arm, respectively. 
Here we refer to the spiral model in Taylor & Cordes 
(1993; Fig. 1 therein) and adopt the circular rotation 
model with Rq — 8.0 kpc (Reid 1993) and the constant 
circular velocity Vb = 220 km s _1 . The tangent point 
in this direction is at ~ 7.0 kpc at the LSR velocity 
~ 103 km s _1 . The HI observation reveals prominent 



absorption features at 68, 81, and 95 km s _1 (e.g., 
LT08). Because these absorption features are seen at 
velocities greater than the systemic velocity, 54 km s _1 , 
of the MC associated with Kes 75, the SNR is located 
at the far side of the Sagittarius arm. This location 
is consistent with the lower limit, ~ 5.1 kpc, placed 
on the distance to this remnant by LT08 based on HI 
absorption. Thus, we determine the distance to be 
10.6 kpc. As a consequence, the 60-110 km s _1 MCs 
are all located in front of the 54 km s _1 MC/Kes 75 
association. 

Our distance estimate is larger than the previous esti- 
mate of the upper limit, ~ 7.5 kpc, set by LT08, which 
suggests no HI absorption by the CO-line emitting cloud 
around 102 km s _1 in the direction of Kes 75. The HI 
spectra in LT08 were extracted from two small regions 
which are bright in the radio continuum. To inspect the 
absorption in the larger radio-bright field, we made two- 
dimensional analysis of HI spectra along the Galactic lat- 
itude and longitude cuts near Kes 75, which are 3-pixels 
wide and cover most part of the radio-bright region (as 
shown in Fig. [TTj) . We see in the produced HI PV dia- 
grams that, above 95 km s _1 at which the prominent HI 
absorption appears, there is a clear absorption feature at 
100-103 km s _1 . Because of this feature, the HI absorp- 
tion at ~ 102 km s _1 , although not strong, cannot be 
ignored, and this gas (at a distance of either 6.4 kpc or 
7.5 kpc) is thus allowed to be located in the foreground 
of the SNR. 

4.2. SNR physics of Kes 75 

The discovery of the MC-SNR association and the up- 
dated distance to Kes 75 are helpful in understanding 
the remnant's physical properties. 

The SNR has a semicircular morphology in Chandra 
X-ray, 24/Ltm mid-IR, and 1.4 GHz radio wavelengths. 
It might naturally have been suspected that the SNR 
interacts from the southern side, and slowed by, the huge 
MC. A similar scenario was suggested for SNR CTB 109, 
in which the SNR interacts with the western massive 
MC (Tatematsu et al. 1987; Wang et al. 1992; Rho & 
Petre 1 997; Sasaki et al. 2004). In this case, for the 
young (Ej3T3j) SNR Kes75, because of the drastic shock 
deceleration and hence the magnetic field compression 
and amplification, the radio emission might have been 
expected to be strong along the interface between the 
SNR and the MC. Also, because of the conversion of 
kinetic energy to the thermal energy, the X-rays might 
have been expected to be bright near the interface. In 
the CTB 109, the SNR is so evolved (~lx 10 4 yr) that 
the X-ray and radio emission have become faint along the 
border. The enhanced emissions are actually seen along 
the western rim of the middle age (<~ 4 kyr) SNR 3C391, 
where the blast wave has been suggested to hit a dense 
cloud plane (Chen et al. 2004). In Kes 75, however, the 
enhancement in either the radio or the X-ray emission 
does not appear along the northern edge of the apparent 
semicircle (if this edge were supposed as the interface). 

An alternative possibility is that the SNR is located 
behind the MC and the X-ray of the northern half suf- 
fers severe extinction. In the latter case, however, since 
the effect of extinction by the cloud is negligible at the 
1.4 GHz, the northern radio emission should be intrin- 
sically faint. Therefore the SNR might propagate into 
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Fig. 11. — The 1420 MHz continuum image (left) and the position- velocity diagrams of HI emission along the Galactic latitude (upper- 
right) and longitude (lower-right) cuts. In the two PV dragrams, HI emission is shown with contours at 20, 50, 70, 75, 82, 85, 88, 92, 98, 
100, 110, and 120 and 10, 20, 30, 50, 70, 80, 84, 87, 90, 93, 96, 106, and 116 K, respectively, with the dark greyscale levels representing low 
brightness temperature. The dashed lines mark the spatial extent of the SNR. The rms noise of the PV diagrams is 1.15 K. 



a northern low-density region behind the MC, and this 
scenario thus explain why the remnant is also faint in the 
north in X-ray and mid-IR. It appears that faint 24/im 
mid-IR diffuse emission can marginally be discerned in 
the northern region (see FigJ12|l , but of course there may 
be confusion between the 24/im emission from the SNR 
itself and that from the larger molecular cloud region. In 
addition, a short, broken structure seems to be present 
in the 89 GHz radio image (Fig. 3 in Bock & Gaensler 
2005), coincident with the northern patch of the bright 
12 CO blue-wing emission (i.e., a weak CO-arc emission 
probably there, §3.2| . namely, related to the northern 
maximum broadening in the CO PV maps (along the 
thick dashed line in Fig [51 i )3.2p . If this radio structure 
is a part of broken shell on the northern rim of the SNR, 
then this part of shell and the southern partial shell (in 
X-ray, IR, and radio) delineate a roughly circular mor- 
phology of the SNR. 

In the south, on the other hand, the partial SNR 
shell seen in these wavelengths has been shown to be 
coincident with the molecular shell in the blue wing of 
the 54 km s _1 CO emission (Fig|7]). The strong radio 
continuum emission may be generated from the SNR 
shock that is slowed by dense clouds (Frail & Mitchell 
1998) or the blast wave that propagates in the intercloud 
material (Blandford & Cowie 1982). The molecular arc 
moving at a velocity v m of order 10 km s _1 (implied 
by the blue-shifted line broadening which reflects the 
velocity component along the LOS), as a possibility, 
may represent the molecular gas that is perturbed by 
the slow transmitted cloud shock after it is struck by 
the blast wave. The position of the southern molecular 
shell may also indicate that some dense molecular 
gas may have been engulfed by the blast wave. The 
thermal X-ray emission may arise from the heated and 
evaporated gas (together with the ejecta) behind the 
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Fig. 12. — The grey-scale Spitzer 24/^m mid-IR image in loga- 
rithmic scales. 



shocked molecular gas, surrounding the engulfed clumps, 
or just behind the blast wave. In this scenario, there can 
be a crude pressure balance between the cloud shock 
and the X-ray emitting hot gas (Zel'dovich & Raizer 
1967; McKee & Cowie 1975): \An m m H v 2 m ~ 2.3n H kT x , 
where n m is the number density of the hydrogen atoms 
ahead of the cloud shock, and kT x the temperature 
of the X-ray emitting gas (see Table 3; here the 
vnei+powerlaw model is used). Thereby n m ~ 4 x 
10 3 (n ff /3cm- 3 )(A:T :l; /0.8keV)(u m /10kms~ 1 )- 2 cm- 3 . 
This is over an order of magnitude higher than 
the mean density of the 54 km s" 1 main-body MC 
[n(H 2 ) ~ 60^6 cm" 3 , %T2] . The 24/xm IR emission 
may come from the dust grains heated by the hot gas 
and even possibly from the shocked molecular gas (e.g., 
OH and H 2 0). 

In the above scenario, the dense molecular shell 
is pre-existing before the SNR shock arrives. The 
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molecular shell could not be the material swept up by 
the supernova blast wave all the way to this radius 
(r ~ 1.8' ~ 6dio.6 pc). We adopt the observed total 
mass ~ 2.7 x 10 3 d 2 6-^0 an( i the Gaussian-subtracted 
mass ~ 0.9 x 10 3 d 2 6 M Q (Table 2) as the upper and 
lower limits of the mass of the southern molecular shell, 
respectively. If it were swept up by the blast wave, the 
original mean density of the molecular gas (counting a 
quarter of sphere of radius 1.8') before it is swept up 
is no(H 2 ) ~ 60— 180tff 6 cm~ 3 (the lower limit of the 
density is similar to the mean cloud density, tj3.2[) . Based 
on the Sedov (1959) evolutionary law, if the blast-wave 
expansion velocity v ~ 3.7 x 10 3 kms~ is adopted 
from the FWHM of the Si line (Helfand et al. 2003), 
the age of the remnant would be ~ OAr/v ~ 640 yr, 
seemingly close to the pulsar's spin-down time scale 
723 yr (Gotthclf et al. 2000). However, the ex- 
plosion energy E = (25/4^) 2.8n (H 2 )m ff u 2 r 3 ~ 
1.3 x 10 54 [n (H 2 )/100cm- 3 ](w/3.7 x 

10 3 kms" 1 ) 2 (r/6pc) 3 ergs (where £ = 2.026), un- 
reasonably higher than the canonical number 10 51 ergs. 
On the other hand, the amount of hot gas in the 
southern clumps (contributing most to the thermal 
X-ray emission from the remnant), ~ 5M Q f£ !3.3j) . 
is over two order of magnitudes smaller than that of 
the southern molecular shell. This would imply that 
massive swept-up molecular gas has quickly reformed 
after they were dissociated and ionized by the blast 
shock (at a velocity higher than ~ 50 km s _1 , Draine 
& McKee 1993). It seems impossible, however, because 
the timescale of H2 gas formation, > 10 5 4 yr (Koyama 
& Inutsuka 2000), is much higher than the spin-down 
time (i.e., the remnant's age). 

We suggest that the pre-existing molecular shell may 
be the cooled, condensed, and clumpy material swept up 
by the stellar wind of the progenitor, which was in the 
vicinity of the large MC, in supplement to the above sce- 
nario of shock hitting the dense molecular gas. The pres- 
ence of the central pulsar evidences that the Kes 75 SNR 
is the offspring of a core-collapse supernova explosion of 
a massive star. Massive stars may create a cavity with 
their energetic stellar winds and ionizing radiation before 
they explode, and the cavity wall may be an imprint of 
the wind modification left on the nearby environment. 
An excellent existing example of the molecular wind- 
bubble shells may be that recently discovered coincident 
with the ring nebula G79.29+0.46, which surrounds a lu- 
minous blue variable star (Rizzo et al. 2008), although 
it has a smaller radius (~ 1 pc) and is at a higher ki- 
netic temperature. As a matter of fact, in the case of 
Kes 75, a cavity is indeed seen surrounded by a molecu- 
lar shell ( §3.2) Fig [7]). The scenario of the molecular shell 
of wind bubble also agrees with the suggestion of a strong 
Wolf-Rayet wind prior to the type Ib/c supernova event 
of Kes 75 (Chevalier 2005). In Kes 75, there seems not 
to be a red-shifted part of shell that expands backward 
and thus the wind-bubble may deviate from the spher- 
ical symmetry in the back side. To catch the essence, 
some parameters are estimated according to canonical 
bubble relation (Castor et al. 1975; Weaver et al. 1977), 
as a crude approximation, as below. The stellar wind 
bubble of the Kes 75 SNR progenitor has an age ~5x 
10 5 (r/6pc) 5 / 3 L^ 6 1/3 [n(H 2 )/100cm- 3 ] 1 / 3 yr, where L 36 is 



the mechanical luminosity of the stellar wind in units 
of 10 36 ergss~ 1 . The bubble shell expands at a veloc- 
ity of - 7(r/6pc)- 2 / 3 J L3/ 3 [n(H 2 )/100cm- 3 ]~ 1 / 3 km s _1 , 
a component of which projected to the LOS may be sim- 
ilar to the blueshift in the broadened line profile of the 
disturbed molecular gas. In addition, the wind-cavity 
wall scenario could not only naturally explain why the 
pre-shock molecular gas hitted by the SNR shock is over 
an order of magnitude denser than the main-body MC, 
but also avoid the above very high explosion energy by 
virtue of the low density (e.g., of order 0.1 cm~ 3 for hy- 
drogen atoms) in the cavity. 

Similar scenario has been suggested to explain the 
physical properties of SNR Kes69, which bears resem- 
blance with Kes 75 in that the molecular arcs are well 
coincident with the SNR shells and move at a velocity 
of order 10 km s _1 (Paper I). In addition to Kes 69, a 
series of other SNRs have been suggested to be born in 
pre-existing stellar-wind cavities, with the shocks hitting 
the cavity walls, such as N132D (Hughes 1987; Chen et 
al. 2003), W49B (Keohane et al. 2007), and Kes 27 (Chen 
et al. 2008). 

In combination of the blue-shifted broadening of the 
CO line profile, the SNR is suggested to be located prox- 
imately behind the 54 km s^ 1 MC, and the molecular 
shell (in the broadened blue wing) was formed mainly on 
the near side close to the MC (i.e. toward the observer), 
as schetched in Figure [131 In the north, the remnant 
material expands into a tenuous region except for the 
side close to the MC, and the northern part of the SNR, 
together with the central PWN, is obscured by the MC 
complexes. 




shell 

Fig. 13. — Diagram showing the positional relation of SNR Kes 75 
with the second and third MC complex components. 

The extra obscuration in the northern area by the 
MCs could essentially account for the difference of Nu 
(~ 1 x 10 22 cm" 2 ) between the northern half of the 
remnant (including the central PWN) and the southern 
half (in the case in which vnei+powerlaw model is ap- 
plied), as derived from the X-ray spectral analysis f ij3.3p . 
This explanation is fairly consistent with the difference of 
the molecular column density ~lx 10 21 cm -2 between 
the northern and southern areas for the second (40- 
60 km s" 1 ) molecular component and ~ 5 x 10 21 cm~ 2 for 
the third (60-110 km s _1 ) component (which is located 
in the foreground, see §4.1|1 . 

4.3. The High-energy Emission of PSR J1846-0258 
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McBride et al. (2008) find that the most unusual fea- 
ture of PSR J1846-0258 is its large efficiency in con- 
verting spin-down power (E) into X- and 7-ray lumi- 
nosity (Xo.5-iokeV and L 20 -iookoV, respectively) of the 
PWN. In their work, the SNR distance of 19 kpc was 
used. If 10.6 kpc is adopted as the distance to the 
PSR, the 0.5-10 keV flux of the Kes 75 PWN 4.0 x 
10~ n ergs cm~ 2 s _1 (Hclfand et al. 2003) corresponds to 
a luminosity Ln.5-iokcV ~ 5.4 x 10 35 crgss _1 and the 
X-ray efficiency is Lo.5-iokcV / E « 6.4%. (The 7-ray effi- 
ciency is corrected to T^o-iookcv/-^ ~ 4.8%.) The X-ray 
efficiency is lower than that (9.1%) of PSR B0540-69 
(Kaaret et al. 2001) and similar to that (5.1%) of the 
Crab [here £o.2-iokcV ~ 3.3 x 10 37 ergs s -1 (Kaaret et 
al. 2001) is converted to £o.5-iokcV ~ 2.4 x 10 37 ergs s _1 
for the Crab Nebula with the X-ray photon index 2.1 
(Willingalc et al. 2001)], and thus the efficiency seems to 
be not peculiar among the young pulsars. 

The High Energy Stereoscopic System (H.E.S.S.) col- 
laboration (Djannati-Atai et al. 2007) discovers the very 
high energy (VHE) emission from the sky region of 
SNR Kes 75 and regards it as a point-like source with 
the point spread function (PSF) 5', which is in a quite 
good agreement with the position of the PSR J1846-0258. 
It is noted that L x /Lj ~ 10 for the case of Kes 75, only 
one-third of SNR G21. 5-0.9 and one-twelfth of the Crab 
Nebula. This seems to imply a 7-ray excess in Kes 75 
comparing to other plerions. Since the H.E.S.S. VHE 
observation cannot resolved the PWN from other struc- 
tures of Kes 75, whose apparent size 4' is smaller than the 
PSF, the 7 excess cannot be excluded to arise from other 
sites than the PWN. The discovery of the association of 
Kes 75 with the MC provides a possible interpretation 
for this excess, because a part of the 7-rays may come 
from the decay of neutral pions produced by proton in- 
teractions between the SNR shock and the MC. 

5. CONCLUSIONS 

Millimeter observations of the molecular gas toward 
SNR Kes 75 have been made in SRAO and PMOD. By 
the spectroscopic analysis of the 12 CO (J=l-0) emission, 
in combination with the X-ray, 24/xm mid-IR, and radio 
continuum data, we conclude the main results as follows. 

1. There are several CO emission peaks in a broad 
LSR velocity range from 3 to 110 km s . The 
~ 83-96 km s _1 MCs overlap a large north-western 
region of the remnant and is inferred to be located 
in the foreground of the SNR. The 45-58 km s" 1 
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molecular gas shows a blue-shifted broadening in 
the 12 CO line profile and the intensity centroid of 
the gas sits in the northern area of the SNR. The 
position-velocity map displays maximums of line 
broadening near the edge of the remnant. 

2. In the integrated intensity map of the broadened 
blue wing at 45-51 km s _1 , a cavity surrounded by 
a molecular shell is unveiled, which is coincident 
with the SNR, and the shell in the south follows 
the bright partial shell seen in X-rays, mid-IR, and 
radio continuum. 

3. All the observational features of the broadened blue 
wing provide effective evidences that Kes 75 is as- 
sociated with the adjacent 10 4 M Q 54 km s _1 MC. 

4. The establishment of the association between 
Kes 75 and the 54 km s _1 MC results in a determi- 
nation of the kinematic distance at ~ 10.6 kpc to 
the remnant, which agrees with a location at the 
far side of the Sagittarius arm. 

5. The X-ray, mid-IR, and radio shell- like morphology 
in the south can be accounted for by the SNR shock 
striking a pre-existing molecular shell. This shell 
may represent the debris of the cooled, condensed, 
clumpy shell of the progenitor's wind-bubble prox- 
imately behind the 54 km s _1 cloud. 

6. The discovery of the association of Kes 75 with 
MC can provide a qualitative explanation for the 
7-ray excess of the remnant. 
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